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Abstract 

The m-bonacci word is a generalization of the Fibonacci word to the m-letter alphabet 
A = {0, . . . , m — 1}. It is the unique fixed point of the Pisot-type substitution (p m : — s> 
■ 01, 1 — > 02, . . . , (m — 2) — s> 0(m — 1), and (m — 1) — ¥ 0. A result of Adamczewski implies 

the existence of constants c' m ' such that the m-bonacci word is c^'-balanced, i.e., numbers 
i^h , of letter a occurring in two factors of the same length differ at most by c' m ' for any letter 

a £ A. The constants have been already determined for m = 2 and m = 3. In this 
paper we study the bounds c- m ' for a general m > 2. We show that the m-bonacci word 
is ([fi:mj + 12)-balanced, where k ~ 0.58. For m < 12, we improve the constant c^" 1 - 1 by a 
computer numerical calculation to the value l" " 1 ^ 1 ] . 
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The m-bonacci word is a generalization of the Fibonacci word to the m-letter alphabet A = 
{0, . . . , m — 1}. It is the unique fixed point of the substitution tp — ip m given by the prescription 

-> 01, 1 -> 02, . . . , (m - 2) -> 0(m - 1), and (m - 1) -> 0. (1) 

In particular, for m = 3, we obtain the substitution — > 01, 1 — > 02, 2 — > with the fixed point 

0102010010201010201001020102010010201010201001020100102010102010010201020100- ■ ■ , 

usually called the Tribonacci word. 

The aim of this article is to study a certain combinatorial property of the m-bonacci word for 
a general m. Namely, we examine the balance property, which describes a certain uniformity of 
occurrences of letters in an infinite word. In order to give its rigorous definition, let us precise the 
notation we will use in the sequel. A factor of an infinite word u = U0U1U2 ■ • • G -4 N is any finite 
string in the form w = u.;Ui + i ■ • • u.; +Tl _i for certain i G No, n G N, where \w\ = n is the length 
of the factor w. The language of an infinite word u, denoted by £(u), is the set of all its factors. 
The number of occurencies of a given letter a G A in a factor w is denoted by \w\ a . Clearly, 
J2aeA \ w \ a = \ w \- ^ ne balance property is related to the variability of \w\ a within the meaning of 
the following definition. 
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Definition 1. Let c be a positive integer. An infinite word u £ is said to be c-balanced if 

Ma - Ma < C 

for all factors w,v £ £(u) of the same length and for each letter a £ A. 

The notion of a 1-balanced word (originally referred to as "balanced word" ) has been used by 
Morse and Hedlund already in 1940 [5] for a characterization of Sturmian sequences. Since the 
Fibonacci word (in our notation 2-bonacci word) is Sturmian, it is 1-balanced. 

It was expected and announced in several papers since 2000 that the Tribonacci word is 2- 
balanced [SJ HI [T3] . This statement has been proved in 2009 (in two different ways) by Richomme, 
Saari and Zamoboni [TT]. As for a general m > 2, in 2009 Glen and Justin [7] mentioned "the 
fc-bonacci word is (k— l)-balanced" , but to the best of our knowledge, no proof of this proposition 
has ever been published. 

The m-bonacci words belong to a broad class called Arnoux-Rauzy words. In the last ten 
years, balance properties of Arnoux-Rauzy words have been intensively studied. For the most 
recent results and a nice overview see [3]. 

The works of Adamczewski on discrepancy and balance properties of fixed points of primitive 
substitutions [UH] imply the existence of finite constants c^ m ^ such that the m-bonacci word is c^ m '- 
balanced. Namely, Adamczewski proved that if all eigenvalues of the matrix of substitution except 
the dominant one are of modulus less than 1, then the fixed point of the primitive substitution 
is c-balanced for some c. It is well known (and explicitly shown in our text as well) that the 
substitution defined by ([1]) satisfies the Adamczewski condition. 

In the present article, we approach the problem of determining c^ m ' by refining the matrix 
method used by Adamczewski in [TJ [2] (and also by Richomme, Saari, Zamboni in [11] in their 
Proof 2). Small values of m can be treated numerically. We show that 

• the 4-bonacci word and the 5-bonacci word are 3-balanced but not 2-balanced; 

• for m — 6, 7, . . . , 12 the m-bonacci word is f 12 ^] -balanced, Theorem 13. II 

The approach works for a general m as well. We prove the following theorem. 
Theorem. (Theorem 16. II ) The m-bonacci word is c^ -balanced with 

c ( m ) = L«mJ +12, 
where k = 2 -r^—, — ^xrflh m 0.58. 

7r JO (5— 4 cos x) ln(5 — 4 cos x) 

Our results confirm the bound c = m — 1 proposed by Glen and Justin for all m < 12 and 
m > 29. Moreover, it turns out that the formerly proposed bound c = m — 1 is far from being 
optimal except for a few small values of m. 

Our article is organized as follows: Section [5] explains relationship between balance and dis- 
crepancy and gives a formula estimating the balance constant using spectrum of the matrix M of 
substitution ((T|). In Section [31 we present results obtained by computer evaluation of this formula. 
In Section |4l we show that for estimating the balance constant c we can concentrate on the letter 
only. Sections [S] and [5] are devoted to the proof of the main theorem. Our proof requires very 
detailed information about spectrum of the matrix M; in Appendix we use standard methods of 
calculus to describe this spectrum. 



2 Balance property and discrepancy 

This section describes the main idea that will be later applied to find for any letter a € {0, . . . , m — 
1} upper bound on the letter balance constant 

c a := max{|w| a - |u| a : v, w £ £(u) and |io| = \v\} . 

The derivation of these bounds uses the following two ingredients. 
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• the m-bonacci sequence defined recursively 

To =Ti = ... = T m -2 = 0, T m _i = 1 



and 



T n — T n _i + T„_2 + • • • + T n - m (2) 



for any n > m: 
• zeros (3 = /3q > 1, /3±, . . . , /3 m -i of the polynomial 

p(x)=x m -x m - 1 -...-x-l. 

It is well known that p(x) is an irreducible polynomial, its root (3 belongs to the interval (1, 2), 
and the other roots (conjugates of /3) are all of modulus less than 1. From now on, we order the 
roots Pi, . . . , ftm-i according to their arguments, i.e., 

< argGSi) < arg(/3 2 ) < • • < arg(/3 m _i) < 2ir . (3) 

The m-bonacci word is a fixed point of a primitive substitution. Therefore, density fi a of any 
letter a e A is well defined and positive, i.e., 

Mq = hm M!ik >0 , 

where u[n] the prefix of u of length n. We refer to [5], where the problem of letter densities is 
studied in detail. 

The value /x can be interpreted in the way that the "expected" number of letters a in the 
prefix u[n] is [i a n. A simple consequence of the definition of \i a is the following observation. 

Observation 1. For any e > and for any positive integer N , there exist factors v and w in 
C(u) such that 

\v\ = \w\=N, \w\ a >HaN-e and \v\ a < ^ a N + e. 

Proof. Assume that there exist e > and N > 1 such that for any factor w of length N, the 
inequality |iu| < fi a N — e holds. It means that for the prefix of u of length n — kN, we obtain 
|u[n]| = \u[kN]\ a < (fi a N-e)k. This implies /i a = lim = lim H^Ik < _ ^ 

which is a contradiction. The proof of existence of v is analogous. □ 

The difference between the expected and actual number of letters a defines the discrepancy 
function D a : N — > R; 

D a (n) = \u[n]\ a - Mq" 

for any n£N. 

Lemma 2.1. For any letter a, denote 

A a := supD a (n) - inf D a (n). 

Then A a < c a < 2A a . 

Proof. Let w,v e ^C(u) be factors of the same length such that c a = \w\ a — \v\ a . We can find 
prefixes W and V of u such that Ww and T^7j are prefixes of u as well. Obviously 

Ma ~ |«|a = \Ww\a - \W\ a - \Vv\ a + \V\ a = D a (\Ww\) - D a (\W\) - D a (\Vv\) + D a {\V\) 
< 2supD a (n) - 2 inf D a (n) = 2A a . 

raGN «e N 
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To deduce the lower bound on c a , let us choose e > 0. There exist prefixes of u, say u[ni] and 
u[rt2], such that D a (ni) > sup„ eN D a (n) — e and D a {n%) < inf„ S N D a {n) + e, or equivalently 

|u[m]| a > /i„ni + sup D a (n) - e, 

|u[n 2 ]|a < Va,n 2 + inf D a (n) + e. 

First suppose that n\ > n 2 and put N := ri\ — n 2 . Denote the suffix of u[ni] of length N by 
W. Then W contains at least /i a N + sup neN D a (n) — inf ng N D a (n) ~ 2e letters a. 

According to Observation [TJ there exists a factor W of length N such that \W\ a < /j, a N + e. 
Hence c a > \W\ a - \W\ a > sup„ eN D a (n) - mi nm D a [n) - 3e = A a - 3e. 

The case n\ < n 2 is analogous. □ 

To find the value A a , we apply the method of Adamczewski used in [TJ [2]. Let us first recall 
the notation used in this method. 

Let M be a matrix of the substitution (UJ). Since entries of M are defined as M a o = |</?(&)| a 
for a, b <E {0, 1, . . . , m — 1}, we have 



M = 



( 1 1 1 

1 

1 

\o 



1 1\ 





1 0/ 



By ty(w) we denote the Parikh vector of the word w £ A* , i.e., = (|m|o, |u>|i, • • • , |w| m -i) T . 

The matrix of a substitution helps effectively calculate the Parikh vector of an image w under tp. 
It is easy to see that 

•$r(tp(w)) = M^(w) foranywe^l*. (4) 

Lemma 2.2. For any prefix u[n] of the m-bonacci word u, there exist £ G N and 5q, Si, . . . , 5g G 
{0, 1} such that 

e 

*(u[n]) = ^4A^(0). (5) 

fc=0 

Moreover, for any choice of £ € {0, 1,2,.. .} and 8q, . . . , 8$ G {0, 1}, there exists a prefix u[n] 
of u such that ^ holds. 

Proof. According to result [6], for any prefix there exist words Eg, ^ e, £^_i, . . . ,E±,Eo (e is the 
empty word) such that 

u[n] = V fi{Ei)<f- 1 {E t - 1 ) ■ ■ ■ ip{E 1 )E Q (6) 

and for any k, the word -Bfe is a proper prefix of ip(a) for some letter a G A. 

For our substitution tp, the only proper prefixes of <p(a) are _Bfe = e and D^ = 0. Since the 
Parikh vector of a concatenation of words is the sum of their Parikh vectors, we have 



*(u[n]) = ^4*(/'(0)), 



k=Q 



where S k = 1 if E h = and <5 fc = if Efc = e. Applying formula flU) to *((^ fe (0)), we get ©. 

In general, not all sequences of Eg, Eg-i, . . . , E\, E correspond to a prefix of u. The relevant 
sequences are described by paths in so called prefix graph of substitution. Nevertheless, since for 



our substitution the equality tp m (0) 
gives a prefix of u. 



^-1(0)^-2(0) . . . ^(0)0 holds, any choice of E l G {e, 0} 

□ 
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Knowledge of the Parikh vector $(u[n]) enables us to compute discrepancy D a (n). To make 
arithmetic manipulation more elegant, Adamczewski denotes row vectors 

hf® = (1,0,. ..,0)-«,(l,l,- •-,!), 
/iW = (0,l,... ) 0)-/ii(l,l,---,l), 



M'"- 1 ) = (0,...,0 ) l)- Mm _ 1 (l ) l,...,l), 

and expresses the discrepancy as the scalar product 

D a (n) = /i (a) *(u[n]). 

Verification of the formula is straightforward. 

Now we can formulate the main tool for estimation of c a . 

Proposition 2.3. For any a g {0, 1, . . . , m — 1} and fceN, denote 

5 (a,fc) = |/(0)| o - M a> fe (0)| , 

where \x a is the density of the letter a in u. Then 

g(a 7 k) = Tk+m-a-\ - ^ a+1 T k+m 
sup D a (n) - inf D a (n) = V* |g(a, k)\ 

nGN * — ' 



(7) 



(8) 



(9) 



and 



k=0 



Proof. At first, since g(a, k) is nothing but D a (\<p k (Q)\), equation © gives g(a, k) = h^^((p k (0)). 
Using equation (g]), we obtain ^((^(O)) = Af fc *(()), hence 



3 (a,fc) = /i (a) M fc *(0) . 



(10) 



This expression combined with equations ((5]) and ([7]) gives D a (n) — J2k=o$kg(a,k), where 

+oo +00 

Sf. € {0,1}. Clearly, supD a (n) < ^ g(a,,k) and inf D a (n) > ff( a i^)- According 

ii£N fc=0 * ™ eN fc=0 

g(a,k)>0 g(a,k)<0 

to Lemma [521 any choice of <Vs corresponds to a prefix of u[n], and, therefore, the equalities are 
reached in the previous inequalities. To sum up, 



fe=0 



sup D a (n) - mi D a (n) = V g(a,k) - V g(a, k) = V" \g(a, k)\. 
„ eN n&> ITo 

g(a,k)>0 

In order to prove equation ©, let us observe that 



k=0 
g(a,k)<0 



( Tn \ 






A 


T n -x 






-2 




= A/ 










\T n - 


-m / 



Since (T m _i,T m _ 2 , . . . , T ) = (1, 0, 0, . . . , 0) = (*(0)) T , we get using (HHD 

3 (a, fc) = /^Af fc *(0) = (T m+fc _ 1) T rn+fe _ 2 , . . . ,T fe ) T 



(11) 



5 



It is readily seen that the vector fl = (/3 _1 , /3~ 2 , . . . , /3 _m ) T is an eigenvector of M correspond- 
ing to the dominant eigenvalue 0. Moreover, sum of components of fl equals 1. It is well known 
that a vector fl with these properties is the vector of letter densities, see [5]. It means that for 
any letter a G {0, 1, . . . , m — 1}, the density of letter a is fi a = f3~ 1 ~ a . If we apply this fact to 
and use the relation we find 



g(a,k) = T m+k - a -i - P~ 



Lrp 

-*- 1 



Corollary 2.4. The balance constants of the m-bonacci word satisfy 

+00 



c a <2j2\g(a,k)\ 



□ 



(12) 



fc=0 

for all a G A. 

Proof. The estimate follows easily from Lemma 12.11 and Proposition 

+00 



a < 2A a = 2 sup D a {n) - inf D a {n) = 2 V \g{a, k)\ 



fc=0 



□ 



Remark 1. To estimate the sum ls( a J ty\> we i«7Z use i/ie explicit formula for elements T n 

of the m-bonacci sequence. The characteristic equation of @ is the polynomial p{x) with zeros 
/3 = @o, 0%, . . . , m —i- Hence there exist constants ao, a%, . . . , a m _i G C such that 

T n = oo(3q + aiP" + . . . + a m —i0 m _ 1 . 

The constants do, a±, . . . , a m -i depend on the initial values T$, T\, . . . , T m _i only. A standard 

calculation provides T n — yZT-n 1 tt% \ 07 ', where p' denotes the derivative of the characteristic 

— 4 3— u p \pj) j 

polynomial p. 

Using (J5J), we can conclude with 



30+1 



P'(0j 



-6 k 

\' J 



k-\-m 



(13) 



3 Numerical upper bounds on balance constant 

According to Corollarv l2.41 the letter balance constants of the m-bonacci word u can be estimated 
by the formula 



c a < 



+ 00 



2$>(a,A:)l 



k=0 



for any letter a G {0, 1, . . . , m — 1} and for all m > 2. 

In this section we estimate the expressions 2 J^t^o \9( a i ^) I usm g a computer calculation. The 
calculations are very time-consuming for m above 10, therefore, we confine ourselves to m < 12. 

The calculation is based on the following strategy. We sum up the first n members of 
(l<7(a,fc)l)/S an d estimate the rest of them; 



+ 00 



Iff (a, fc)| < \g(a, k)\ + E, where E satisfies E > \g(a, k)\. 



fc=0 



fc=0 



k—n 
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Table 1: 4-bonacci - g(a,k) with quadruple of integer coefficients in linear combination of 
g(a, 0), . . . , g(a, 3) and its signum. 





IC of (5(o,£;))ft=o 


a = a = 1 a = 2 a = 3 


g(a,0) 


(1,0,0,0) 


+ 




(U, 1, U, UJ 


1 




(0,0,1,0) 


- - + - 


g(a,3) 


(0,0,0,1) 


+ 


ff(M) 


(1,1,1,1) 


+ - - - 


g(a,5) 


(1,2,2,2) 


+ 


5(0,6) 


(2,3,4,4) 


- - + - 


5(a,7) 


(4,6,7,8) 


+ 


5(0,8) 


(8,12,14,15) 


+ + - 


5(0, 9) 


(15,23,27,29) 


+ + 


g(a, 10) 


(29,44,52,56) 


- - + + 


o(o,ll) 


(56,85,100,108) 


+ - - + 


5(o,12) 


(108,164,193,208) 


+ + - - 



Formula (flU)) provides setting 

i-i&r 

To conclude, we have to find an n big enough to satisfy 

(14) 

Since we always compute on machines working in a finite precision, it is desirable to reduce 
the work with non-integer numbers. Therefore, we make use of the fact that, for a fixed letter a 
and the alphabet cardinality m, the sequence of numbers g(a, k) satisfies the m-bonacci recurrence 
relation 

5(0, n + m) = g(a, n + m — 1) + . . . + g{a, n) , 

which follows from Proposition 12. 31 

Let us demonstrate the method on the 4-bonacci word. The first step is calculating sgn g(a, k) 

from ([9]) for all k £ {0, . . . , m — 1} (illustrated in Table [lj. Then we express X)fc=o \9( a ' &)l as an 

/ 9(0,0) \ 

integer combination (IC) of I : J , which can be rewritten in the form p + „^ +1 for some 

\ </(a,m-l) / 

p, q G Z (this follows from Proposition I2.3P and then evaluated 1 (see Table [5]) . The final step is 
verification of the equality (fill . 

To make our procedure reliable with respect to possible rounding errors, we replace the esti- 
mated error E a ^ m by a constant E > E a , m . If (IT4l) holds, it is equal to the desired upper bound of 
c a (but it may not be optimal). In the opposite case, we must increase n and repeat the procedure. 

Our results obtained for m £ {2, . . . , 12} are summarized in Table [3] 

To find lower bounds on the constant c, one needs to find two factors v,w of the m-bonacci 
word that are of the same length with |w| a — |i>| a big enough. Computer searching in the set of all 
factors is very time-consuming. Nevertheless, for any given m > 4 and any a £ {1, . . . , m — 1}, a 
modification of the abelian co-decomposition method 12 allowed us to find a pair of factors v, w 
of the m-bonacci word such that \v\ — \w\ and \v\ a — \w\ a = 3. For instance, if m = 4, the words 

t, = V 2 (0)^ 9 (0)<A0)^(0), 

1 The calculation must be performed in an environment working in enough precision, e.g., Wolfram Mathematica. 



E a ,n := |/3 2 | n £ 



r 1 



1 



. k=Q 



]{a,k)\ 




g{a,k)\ + E atJl 
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Table 2: 4-bonacci - Estimates of J2k=o If ( a ' ^)l an( ^ the resulting upper bound on c a . 





a = 


a= 1 


a = 2 


a = 3 


12 

^2\g(a,k)\ as 

fe=0 

IC 


/123\ 

1 215 1 

V 232 / 


(1) 


/-133\ 

1 -233 1 
V -254 / 


(=1) 

V -86/ 


12 

k=0 

symbolic 


1664 - ^ 


286 


^-487 


^-86 


12 

^2\g(a,k)\ nu- 

fe=0 

merical 


1.2778 


1.5157 


1.5611 


1.5776 


E a ,13 


0.20054 


0.22213 


0.25916 


0.31056 


12 

£>(a,fc)|+£ 

fe=0 


1.49844 


1.76006 


1.84618 


1.91919 


c a upper bound 


2 


3 


3 


3 



Table 3: Upper estimates of c a forme {2, . . . , 12}, a £ {0, . . . , m — 1}. 



m \ a 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


2 


1 


1 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


3 


2 


2 


2 


X 


X 


X 


X 


X 


X 


X 


X 


X 


4 


2 


3 


3 


3 


X 


X 


X 


X 


X 


X 


X 


X 


5 


2 


3 


3 


3 


3 


X 


X 


X 


X 


X 


X 


X 


6 


3 


3 


4 


4 


4 


4 


X 


X 


X 


X 


X 


X 


7 


3 


4 


4 


4 


4 


4 


4 


X 


X 


X 


X 


X 


8 


3 


4 


4 


4 


4 


4 


4 


4 


X 


X 


X 


X 


9 


3 


4 


5 


5 


5 


5 


5 


5 


5 


X 


X 


X 


10 


3 


5 


5 


5 


5 


5 


5 


5 


5 


5 


X 


X 


11 


4 


5 


5 


6 


6 


6 


6 


6 


6 


6 


6 


X 


12 


4 


5 


6 


6 


6 


6 


6 


6 


6 


6 


6 


6 
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w = (^ 9 (0V 8 (0V 5 ((V(0)) V 11 (00)^ 10 (0)^(0)^(0)^(0)^ 3 (0)^(0)0 
are factors of u such that \v\ — \w\ — 3305, — |w|i = 3. Similarly, if m = 5, the words 

v= V 4 (ov n (ov 6 (oV(o), 

w = (^(O^m^iO^iO))- 1 ^ 13 (00)^ 12 (0)^(0)^(0)^(0)^(0)^(0)^ 2 (0)0 

are factors of u such that \v\ — \w\ = 15481, 1 1 1 — |u>|i = 3. 
Therefore, we can conclude with the following theorem. 

Theorem 3.1. For m G {4, 5}, the m-bonacci word is c-balanced with c = 3 and this bound cannot 
be improved. 

For m G {6, . . . , 12}, the m-bonacci word is c-balanced for c = l" "^ 1 ] . 

4 Balance property of letters in the m-bonacci word 

The numerical calculation, performed in Section [3l is convenient only for small values of m. In 
the rest of the paper we develop a technique to estimate the constant c for the balance property 
of the m-bonacci word for a general m. The calculation will be again based on formula (fT2"j) . but 
this time we bring in an improvement. Instead of estimating the sums X}fe==o ls( a i ^)l f° r au letters 
a G A, we show that in case of the m-bonacci word, the balance constants c a for a = 1, 2, . . . , m — 1 
can be estimated by a simple formula in terms of Co providing that c$ is small enough, see the 
following observation. 

Proposition 4.1. Let m > 4. If c < 2 m_1 - 3, then 

C ,<(2-^)co+4(l-l) (15) 

for each j = 1, 2, . . . , m — 1. In particular, the m-bonacci word is c-balanced with c = 2cq + 3. 

With regard to this proposition, it will be sufficient to estimate ls( a '^)l ano - use f° r ~ 

inula (TT21 just once, for a — 0. All the remaining constants c a for a — 1, . . . , m — 1 can be then 
easily estimated using formula (|15l) . 

Before we prove Proposition 14. 11 we derive two simple observations. 

Observation 2. For any factor f of u and for each j G {1, . . . ,m — 1}, it holds 

|/|o = \<p>{f)\i and |/| = \<f?{f)\i-i. 

Proof. From the form of the substitution (flj, we see |w|j^i = |^(t/;)|j and |u>| = |<p(w)|o for any 
factor w and letter j — 1, 2, . . . , m — 1. Applying these relations on w = /, w — (f(f), ■ ■ ■ , w = 
L P^~ 1 (f)> we get the formulae in the observation. □ 

Observation 3. If f is a factor of u such that \ f\ < 2 m , then |/|o < ^|/| + 1. 

Proof. The form of the substitution ip implies that 00 is the longest block of zeros occurring in u. 
Further, with exception of this block, the letter is always sandwiched by nonzero letters. It is 
easy to see that the shortest factor w ^ 00, with the prefix 00 and the suffix 00 such that w has no 
other occurrences of 00, is the factor w = 0ip m (0)0. Since \w\ = 2 m + 1, any factor f with |/| < 2 m 
contains at most one block 00. This implies the inequality for |/|o stated in the observation. □ 

The following lemma is the combinatorial core for the proof of Proposition |4Tj 

Lemma 4.2. Let j G {1, . . . , m - 1}. // Cj_i < 2 m - 2, then 

c 3 <c + 2 + ^-. (16) 
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Proof. With respect to the definition of Cj, there exists a pair of factors v and w such that 

H = H and \v\j-\w\j = Cj. (17) 

Without loss of generality, we can assume that v and w is the shortest possible pair satis- 
fying (|17|) . Then v and w are in the form v = and w — £■■■£' for certain £,£' ^= j. 
Moreover, we can assume that jw is a factor of u (otherwise we replace w = Uj • • • u^i^i-i by 
w' = Ui^i' ■ ■ ■ VLi+\ w \-x-i' without violating equations (|17p). 

Because of the form of v, there exists a factor V = OV G C(u) such that v — j^(V'). Clearly, 
v is a suffix of <p>(0V) = ^'(V). 

Let wzj be a factor of u such that \z\j =0 (we extend the factor w to the right up to the next 
letter j). As jwzj £ £(u) by assumption, there exists a factor W such that wzj = (pi (WO). 

Observation [5] implies 

. |y| = i + |V|o = i + W 3 iv% = W(v% = \v\i 

. |W] = \W% - 1 = 1^(^0)1,- - 1 = \wzj\j - 1 = \w\i 

• \v\ = i + = i + \^(y%-i = i + \v\j-! 

. \W\ = \W0\ - 1 = ^{WO^j-l - 1 = \wzj\j-! - 1 
Together, we have deduced 

\V\o-\W\ =c 3 and \V\ - \W\ = \v\j-! - K'_i - \z\j-i + 2. (18) 
We distinguish two cases: 

• Case \V\ < \W\. Let V = Vx be a factor of u such that \V\ = \W\. From the definition 
of c and C0§ we get c > |F| - |W|o > |^|o - \W\ = cj. Thus Cj < cq + 2 + ^ holds 
trivially. 

• Case |V| > |W|. Let W = Wy be a factor of u such that \W\ = \V\. Then c > | V| — |^|o = 
| V|q — |T4^|o - \y\o = c 3 — \y\o due to ([15]), To bound length of y, we apply Equation (fT5)). It 
gives |y| = |V| - |W| = - - + 2 < - |*ujj-_x + 2 < c,-_ x + 2. With 
regard to the assumption Cj-\ < 2 m — 2, we have |t/| < 2 m . Therefore, |y|o < + 1 due 
to Observation [3] To sum up, cq > Cj — (|(cj_i + 2) + l) >Cj — 2— \cj~\. 

□ 

Proof of Proposition ^. 1\ Let us assume that cq < 2 m ~ 1 — 3. We prove equation (fT5|) by induction 
on j. 

I. Let j = 1 . It holds Co < 2™ l ~ 1 — 3 < 2 m — 2 by assumption, therefore, we can use Lemma (|4.2[) . 
It implies ci < Co + 2 + ^, hence indeed Cj < (2 — A) c + 4 (l — jr) . 

II. Let j > 1 and equation (|15p /ioZd /or j — 1. Inequality Co < 2 m_1 — 3 implies 

1 \ M . . .-^ m _x 



c,--i < - — T J c + 4 ^1 - — T j < 2c + 4 < 2(2 m - J - 3) 
It allows us to apply Lemma (|4.2[) . Equation (fl7j)) gives 



c i <co + 2 + ic i _ 1 <c + 2 + i((2-^ r ) Co+4 ( 1 -2^)) = ( 2 -^) ^(^^J 

In particular, (1151) yields Cj < 2co + 4. As c = max{cj : j = 0, 1, . . . , m — 1} and c and Co are 
integers, necessarily c < 2c + 3. □ 
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5 Estimate of £™ \g(0, k) \ 

As anticipated in Section^ the balance constant cq will be obtained using formula 11 21 Therefore, 
we need to estimate the sum J^t^o This is the topic of this section; since we deal with 

the letter a = only, we abbreviate the symbol <?(0, fc) to g(k). 

The sum J^t^o |p(0, k)\ will be estimated by splitting it into two parts, X)fe=o IsCOl an d 
SfcSm Ifl'WI' an< ^ estimating each of them separately. In Sections |5 . 1 1 and HT21 we show that 

2m— 1 p. +oo ^ 

^2 < 4 an d ^ WOI < ^ m + 1 for all m > 4. 

fc=0 fc=2m 

To get these estimates we will exploit bounds on absolute values and arguments of zeros of poly- 
nomials p{x), derived in Appendix IA1 



2m- 1 

5.1 An upper bound on the sum \g{k)\ 

k=0 

At first we express g(k)'s for all k = 0, 1, . . . ,2m — 1 and determine their signs. Recall that 
fio = 1/P, therefore, due to equation @, it holds 

ff(*)=|^(0)| - ~|/(0)| • (19) 

In the sequel we use the following formula to calculate g(k) for all k < 2m — 1. 
Proposition 5.1. It holds 

I = / 2 fe for k = 0,...,m- 1; , 

IV Wl | 2 fc - 2 fc -" 1 - (k - m)2 k -" 1 - 1 for k = m, ... ,2m - 1 . 1 ; 

Proof. The identity /(0) = y> (/ _1 (0)) together with the substitution (JTJ implies 

|/(0)| =2|/- 1 (0)|-|^- 1 (0)| m _ 1 . (21) 

Let us distinguish two cases. 

• Case k < m — 1. It holds /(0) = and |/(0)| m _i = for all k < m — 2, hence 
1/(0)1 = 2|/- 1 (0)| for all k < m - 1. 

• Case fc > m. We prove equation (|20p for k G {m, m + 1, . . . , 2m — 1} by induction on k. 

I. k — m. We have |v9 m_1 (°)U-i = 1, hence |<^ m (0)| = 2|^ m ~ 1 (0)| - 1 = 2 m - 1. Since 
2 m _ i _ 2 ™ _ 2" 1 -™ - ( m - m)2 m ~ rn ~ 1 , the statement holds true for k = m. 

II. k > m + 1. Let |/ _1 (0)| = 2 fc ~ 1 - 2 fc - 1 " m - (fc - 1 - m^" 1 -" 1 " 1 . The identity 
|/ _1 (0)| m _i = \tp m (0)|, valid for every fc > m + 1, allows us to use the formula (I2T1) in 
the form 

|/(0)| =2|/- 1 (0)|-|/- 1 -™(0)|. 
Since fc — 1 — m < m — 1, we can apply the results obtained above fc < m — 1, whence we get 
|/(0)| = 2 (2 fc ~ 1 - 2 k - 1 - m - (fc - 1 - m)2 k - 1 -' m - 1 ) - 2 k - 1 - m = 2 k - 2 fc ~ m - (fc - m)2 fc -™- 1 . 

□ 

To determine signs of g(k)'s defined by (|T9|) . we need a fine estimate on j3 . Let us recall that 
P is the dominant eigenvalue of the matrix of substitution M and thus a zero of its characteristic 
polynomial p(x) — x m — x" 1 ^ 1 — x m ~ 2 — ■ ■ ■ — x — 1 . 



11 



Proposition 5.2. It holds 
<7(0) = l-~>0; 

g(k) = 2*- 1 (l - |) < /or fc = 1, . 

,(m) = 2— (i-|)+i>0; 



,m-l; 



5 (fc) = f 1 - (2 fc - 1 - (fc + 1 - m)2 fe - m - 2 ) + i • 2 fe - m - 1 < for k = m + 1, . . . , 2m - 1 . 

Proof. The formula for g(0) follows immediately from equation (|19l) . 
For every fc > 1, it holds |<£ fc (0)| = \ip k ~ 1 (0)\, hence 

ff (fc) = |/-i(0)|-i. |/(0)|, 

cf. equation (TK)1) . All the formulae for g(fc) listed in Proposition 15.21 then follow easily from 
equation ((20"]) . 

In the rest of the proof we show that g(0) > 0, g(m) > 0, and g(k) < for all k £ {1, . . . , m — 
l}U{m + l, ...,2m- 1}. 

At first, (3 £ (1, 2) immediately implies g(0) > and g(k) < for all k £ {1, . . . ,m — 1}. 

Ac frw Z" — >m cVioll cVir^iTT 4"Vio4- 



As for k — m, we shall show that 

>m-i ( 1 _ _ 

/?7 /3 



l-^+i>0. 



This inequality is equivalent 



to 



2 - /3 < ' 



2m— 1 ' 



which is valid due to (031) from Appendix, because l/2 m 1 > l/(2 m - (m + l)/2) for all m > 2. 
Similarly, if fc > m + 1 , we need to prove that 

1 - | J (2 fe - 1 - (fc + 1 - m)2 k - m - 2 ) + ~ ■ 2 k - ,n - 1 < for k = m + 1, . . . , 2m - 1 ; 



i.e. 



2-B> i for all fc = m + 1, ... ,2m- 1. 

Z 2 

Since fc + 1 — m < m, the validity immediately follows from inequalities (|43[) . □ 
Proposition 5.3. holds 

2m-l /<? \ 1 1 

E = 1 + (j - 1 j t 2 " 1 (^^ - 1) - (m - 1)2™- 2 ] - - (2'- 1 - 1) < 1 + - . (22) 

Proof. Proposition 15.21 implies 

2m— 1 m — 1 2m — 1 

^|#)h S (0)^#)+ 9 (m)- £ 5 (fc). 

fc=0 fc=l fe=m+l 

When we substitute for g{k) from Proposition 15.21 we obtain 

.g(0)+.g(m) = l + 2'"- 1 (l - -0 , 
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E^h-E 2 ^ 1 i 



fc=i 



_^m-l _ 1) 1 



fc=l 

and, in a similar way, we get 

2m- 1 / 9\ 1 

- E -9( fc ) = - )[2 m (2 m - 1 -l)-(m-l)2 m - 2 ]--(2 m - 1 -l) . 

fc=m+l VP/ P 

Summing up these expressions, we get formula (|22[) . 

In the rest of the proof we show that Y^t—o 1 < 1 + 1/4, which is obviously equivalent to 

(2 - (3) [2 m (2™- 1 - 1) - (to - l)2 m ~ 2 ] - 2" 1 - 1 + 1< J , 



and also to 



(2-/3) 



2 m £ 2 m-l _ 1) _ ( m _ 1)2" 1 - 2 + 



)fTt-l 



1 < 



Using inequality (l4*3l . we obtain 

(2-/3) 2 m (2™- 1 - 1) - (to- l)2 m ~ 2 
1 



< 



Om _ m+1 
Z 2 



2 m (2 m_1 



1 



-i m+l 
~ 2^TT 



- 2" 1 - 1 + 1 = 



r 

4_ 


_ 2 m-l + j 












1) 


- (to - l)2 m - 2 + 


1" 

4_ 


- 2 m_ 




1 




m — 1 i 1 i m+l 




m+l 


1 


1 - 


2m+l 




4 t 2m +a -r 4 




2 m+l 


2 ™ + i 




1 m+l 
-L 2 ™+! 






2 


1 - 


m+l 
2 m+l 



< 



□ 



5.2 An upper bound on the sum Yl 

k=2m 

Proposition 5.4. For any k € N we have 



\g(k)\ < 



m— 1 



2(m- 1) 



ETO)-i|-lftf • 



(23) 



Proof. With regard to equation (|42|) from Appendix, 



.T - 1 



... (to + l)x m - 2mx m ~ 1 x m+1 - 2x m + 1 (to + l)x m - 2TOX" 1 " 1 p(x) 

p W = 1 ? T\2 = 

x — 1 (a; — 1J Z x — 

Since p(/3j) — for every eigenvalue of M, we have 

(to + l)/3™ - 2to/3™ - 1 ( m + l)^ - 2to 
P (ft ) = i ; = tt^- 

Therefore, due to ([lj 



- 1 



0,-1 



m — 1 



l l 



g(k) = E 

As g(k) is real, we can write 



Cjk-\-m 



m — 1 



# y (m+l)/3j — 2m /^m 

■■=i VP - ? p/ Pj 



= v g - ft ft - i g fc 

(m +!)#,- - 2to Pj ' 



(ft - 1)# ) , 



(24) 
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and estimate 

m— 1 1 



(m + l) ( s/-2m (/3j ' 



m — 1 1 

< v - 



(to + l)Pj - 2m 
To finish our proof we will deduce for all j = 1, . . . , m — 1, 

/3 - ft < 1 

" 2(to- 1) ' 

TO — 1 — (TO — 1)t^- 



TO) " II ' \$ 



(m + l)pj - 2m 



Since 



P-p, 



(m + 1)j8j - 2to 



1 



2(m- 1) 



TO — (TO + 1) 



(25) 



(26) 



it suffices to prove that 



TO — 1 — (to — l)4p 



to — (m + l)^f 



< 1. 



We have 



to — 1 — (to — 1)4- 



to — (m + l)^f 



to — 1 



^(ft) 



-1 2 



[ m _r^ K(A .^ + [m±l^.)]- 



Lemma IA.1I implies 2 — /3 < < ^h- ; hence 



TO — 1 TO + 1 
^ < " 



Therefore 



to — 1 



In what follows we demonstrate that 

to — 1 



< 



TO + 1 



3(ft) 



TO — 1 



< 



(27) 

(28) 
(29) 

(30) 



Since P £ (1, 2) and \pj\ < 1, we have 



to — 1 , v to + 1 , 
< to - 1 — 5ft (ft) = to ^— RCSj) - 1 



m + 1 to — 1 



It holds $t(Pj) < 1, and the expression — ^-g^ is positive due to equation |2"51 therefore 



m + 1 to — 1 



< -1 + 



TO + 1 TO — 1 



-(to — 1) 



1 1 

P~ 2 



< 0. 



Hence 



< to - 1 — < m jf- »G8j) . 



£ 2 
i.e., (|3T)]) holds true. Equation (f2"T| together with inequalities (f2T)| and (|3ll|) implies 

2 



to — 1 — (m — 1) — 



to — (to + 1) 



ft 



< 



[2m - (to + l)g(^)r + [(to + l)S(ft)] 2 
[2m - (m + l)3?(/3,)] 2 + [(to + l)3(ft )] 2 



1 . 



□ 



14 



Corollary 5.5. 



+00 m-l 



Proof. Using <|25)l . we can estimate 

+00 1 m — 1 +00 

i V 

2(m - 1) ^ 



(31) 



L\J lit — _l_ ^JV-I I I L — _l_ 



fc=2m 



fc=2m 



Finally, we use Observation 0] to rewrite |/3j| 2 " 1 = 1/|2 



□ 



At this stage we apply the information on \/3j\ for j — 1, . . . , m — 1, derived in Lemma I A. 2 1 
Proposition 5.6. If /10/ds 



+00 ^ / 2 m m_1 

^ lff(fc)l < 2(m-l) , l-21n3 J3¥f E ' 
fe=2m v ; 3m \ m 7 = 1 



1 — COS 7j 



771 — 1 

E 



cos 7j 



-J (5 — 4cos7j) ln(5 — 4COS7J-) ^— J 5 — 4cos7j 

(32) 

Proof. We will estimate summands from inequality (|3Tj) . In the notation /3j = Bje 1Jj , we have 
|SR(/3,-) - 1| _ 1 - cos 7j _ 1 - cosTj- 



1 -B, 



1 - £L 



cos 7j 



thus equation (|44|) from Appendix implies 



1 — cos 7, 1 — cos 7, 2m 

— — + COS 7» < —7- - • 5—5- + COS 7, 

1-Sj J ~ ln(5 - 4cos7j) 1 - Ml ' 3 



(33) 



Concerning the term 1/|2 — /3j| 2 , it holds 



|2 - /3j| 2 4 - 4B 3 - cos 7j + Bf 5 - 4cos7j + 4(1 - Bj) cosTj - 2(1 - Bj) + (1 - B 3 ) 2 



< 



5-4COS7J 1 _ (1 _ Bj) 2 ' icosli 



5—4 cos 7j 



It is easy to see that 2 7 cos 7 < ~, therefore, it suffices to estimate 1 — Bj from above. Since 

J 5 — 4 cos 7 — 3' 1 J 



and 



it holds B 3 > ^ 



V3 = e 1 ^ < 



> 



1 1 



4 - ABj cos 7j + B 2 2 V3 



1 \ >° 3 



m 1 

In 3 1 



In 3 



m 1 

In 3 1 



Hence 1 — Bj < ^ for all j = 1, . . . , m — 1. Consequently, 



< 



|2-/3j| 2 " 5-4cos7j 1 _ 2 . ins • 



3 m 



Inequality f|3T|) combined with estimates (|33l) and f|34)l leads to formula ([32]) . 



(34) 



□ 



The following lemma is an essential component of our calculation. It uses the information on 
7j obtained in Lemma I A. 3 1 
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Lemma 5.7. It holds 

TO— 1 -, /.2-7T i 1 1 / i 

El — cos 7,- to / 1 — cos x , 1 m-1 i / 1 

._ (5 — 4cos7j)ln(5 — 4cos7j) ~~ 2ir J (5 — 4cos7j) In (5 — 4cosx) 6 m 16 \ 36 

(35) 

Proof. Let us denote 

N 1 - cos X 

fix) 



ln(5 — 4 cos a;) ' 

then 



^ l-cos 7i TO ^ 2tt 

^ ln(5-4cos7 J ) = ^ ^ ™ 



The estimate implies 7^ e (j - |) , ^ (j + |)). Therefore, the sum fl3H) is a Riemann 
sum of the function / with respect to the tagged partition 

7T 7T 27T / l\ 

— = xq < x% < ■ ■ ■ < ^m-i = 2-7T , where a;,- = — \ j + - , 

to to m \ 2 y 

of interval [— , 2ir — — 1 . Let us rewrite the summands of (151)1) using a trivial identity 

2tt r x i r Xj 

— fiij) = / m dx + / (/( 7i ) - /(£)) dx . 



in 



Xj-l 



Since 



/(7i)-/(s)<l*-7i|- max {|/'(x)|} < |z - 7j| • max {|/'(x)|}, 



we have 

— f(7j)< r f(x)dx+ max P \x- l3 \dx. 

Now we apply another identity, valid for any jj G [xj-i, Xj], 

\x — 7j|da;= / (7, — x) dx + / (x — jj)dx= / xdx+ xdx. 

Jxj-i Jjj Jo Jo 

which provides us, using the estimate (l52l) . the inequality 

r^. + -sk r-Sc—sk 7T 2 / 1 

\x — 7j j dx < / xdx + xdx = — - 1 



36 



Hence 



27r r Xj 7T 2 / 1 

-/(7i)< / f(x)dx+ max {|/'(x)|}— 1 + - 

TO Jx _i xe[0,2-7r) TO \ 36 



Consequently. 

m-1 / „2ir 



£/( 7i ) J . m /(x)d, + (m-l) 4o ax { |/'(,)|}^(l + 

J = l \ rn 



36 



Furthermore, it can be checked that /(x) > l/6for all x G (0, 7r/2) U (37r/2, 27r) and limx^o /(&) = 
1/4 > 1/6, hence 

JJ m f(x) dx = JJ f{x) dx - J m f{x) dx - jj ^ fix) dx < JJ fix) dx - ^ • i . 
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Finally, a numerical calculation gives max^g [0,277) {|/'(^)|} < g- To sum up, 

m — 1 

El — COS 7j 
_ x (5 - 4cos7j) ln(5 - 4 cos 7^) 

to ( /" 27r 1 - cos x , 2tt 1 , ,1 i 2 / 1 



-2w\J (5-4cos7j)ln(5-4cos:c) dx to ' 6 + ^ ^8 ' m 2 [} ' 36 
whence we obtain the sought formula (l35|) . □ 
Lemma 5.8. It; /ioWs 

m— 1 r 

^ 5-4cos7, " 6 + 6' ld7j 
j=i 13 

Proof. If we define j m := 27r, we can write 

7J / — 1 772 

ST-\ COS 7j m 2lT COS 7j _^ 

^ 5 — 4 cos 7, 27r ^ m 5 — 4 cos 7,- 

3=1 ' J 3=1 ' J 

The sum 2j=i ^f(lj) f° r / (7) := 5-4COS7 will be calculated in a similar way as in the proof of 
Lemma 15.71 Namely, it is a Riemann sum of the function / with respect to the tagged partition 

7T 7T 2ir ( 1 \ 

— = x < xi < . . . < x m -i < x m — 2n H , where x, ; = — \j + - , 

m to to \ 2/ 

of interval 27r + ^] . Following the steps of the proof of Lemma T5.71 we obtain 

- TO / r' 27r+ ™ IT 2 f 1 \ 7T 2 N 

£/(7i)<£M / /(x)dx + (m-l) max ||/'(x)|} — 1 + - +— max {|/'(x)|} 

^— J 27T yJ^L 2;G[0,2tt) TO^ \ 36/ TO :r£[27r,27r+7r/m) ^ 

< ^ f / 2T+ " /(*) dx + to • max {|/'(x)|}4 f 1 + ^ 
2tt \ Jjl. xe[o,2n) to \ 36, 



With regard to the properties of cos, we find 



s: 



cosx , „ / cosx 

dx = 2 / dx = 2 



5 — 4 cos x ./n 5 — 4 cos x 

Furthermore, 



1 — arctan ( 3 tan — 

4 6 



7T 

3 



-x { |/' ( x)| } = I . V^i^-n < | 

ze[o,27r) 2 (15-7T53) 2 8 



To sum up, 



771— 1 



Ecos7j m { ir 9 7T / 1 \\ to 7r / 1 1 \ 9 to 5 

_ x 5-4COS7J - 2^\3 + 8'm\ 1+ m))^ 1 ^J + 2\ 1+ mj8^ 1< ~6 + 6 



□ 



Proposition 5.9. For all m > A, it holds 

E i-9( fc )i < £ m+1 ' ( 38 ) 

where 



2tt 

k—2m 



A := [ - 1 -di» 0.909. (39) 

7 (5 - 4 cosx) ln(5 - 4 cosx) 
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Proof. Recall that 

+00 ^ -tin m—1 m— 1 

V \a(k)\ < m V l-cos 73 + V C ° S7j 

^ m ;l 2(m- 1) 1 - 2inl I j _ in3 (5 - 4 cos 7,-) In (5 - 4 cos 7, ) ^ 5-4cos7, 

cf. formula ([52)1 . If we estimate the sums using inequalities (|35D and (1571) . we obtain 

+00 

\n(k\\- 

2tt 



]T | ff (fc)|-i. m _i 



fe=2m 



1 1 ( 2m fm 1 m -1 ir / 1 \\ m 5\ A 

< 2(m - 1) ' 1-^1 1 1 — las ^ A ~ 6 + "r^T" 'l6l v 36yJ + 6" + 6) _ 2^™ ~ ' 



A numerical integration gives A ps 0.909 £ (0.9, 0.91). For such value of A, the expression above 
is negative for all m > 4; i.e., 



+2^ A 



□ 



6 Main result 

Theorem 6.1. For every m > 5, £/ie m-bonacci word is c-balanced with 

c = [_«7TlJ + 12, 

wAefC K = f J,, 2 " ^cos^ln^cos.) ^ * °- 58 ' 

Proof. In Propositions 15.31 and 15.91 we showed 

2m— 1 +oo ^ 

^ |. 9 (fc)| < - and ^ < — m + 1 for all m > 4 ; 

k=0 fc=2m 

therefore, 

El^)! < 3 + ^™' ( 4 °) 
where A = [^ x \ c ° s ,? A -Ax w 0.909. 

JO (5— 4 cos x) ln(5— 4 cos a;) 

Having this bound in hand, we can use Corollary 12 .41 to estimate the balance constant of letter 
by 

+2° A 
co<2V|g(0,fc)|<| + -m. 

Since § + — m < 2 m ~ 1 — 3 for any m > 5, the assumption of Proposition l4.il is fulfilled and 
thus the m-bonacci word is c-balanced with 

+°° 2A 
c = 2c + 3 < 3 + 4 > \g(0,k)\ < 12+ — m, 

* — ' 7T 

fc=0 

which proves the theorem. 

□ 
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A On eigenvalues of M 

In this section we examine the eigenvalues of the matrix of substitution. In particular, we esti- 
mate their absolute values and arguments. Such information is essential for estimating the sums 
E^o" 1 and £+r 2m |fl(0,*)| in Section^ 

Let us recall that the eigenvalues of the matrix of substitution M are zeros of its characteristic 
polynomial p(x) — x m — x" 1 ^ 1 — x m ~ 2 — . . . — x — 1. The following observation will make further 
calculations substantially simpler. 



19 



Observation 4. Every zero of the polynomial p{x) is a root of the equation 

x m (2-x) = l. (41) 

Proof. For every x ^ 1, we can write 

p(x) = x m - ^— ^ = _ + . (42) 

In particular, p(f3j) = implies — 2/3™ + 1 = 0, whence f}j is a root of equation (14T1) . □ 

At first we derive a fine estimate on j3, which is needed for calculating the sum X)fc=cT lff(0> 

Lemma A.l. The dominant eigenvalue [3 > 1 of the matrix of substitution M obeys the inequal- 
ities 

2m~~~m < 2 - P < 9m m+1 ■ ( 43 ) 
Z 2 Z 2~~ 

Proof. Observation 2] implies /3 m (2 — /3) = 1, hence /3 < 2. Let us set xo := 2 — (3. Obviously, xq 
is a root of the polynomial 

q(x) = (2 - x) m ■ a? - 1 . 

Since /3 s (1,2), necessarily xo € (0, 1). It holds (x) = (2 — a;) m_1 (2 — a; — mi), therefore, g 
grows in [0, 2/(m+ 1)] and decreases in [2/(m+ 1), 1]. Since q(0) = —1 and q(i) = 0, the root Xo 
belongs to the interval (0,2/(m+ 1)), in which q grows. Consequently, proving inequalities (|43l) 
consists in showing that 

l A ( l 



g^ 2 "- f J <0< ^ v 2». BLfal 

Let us start with the estimate of 2 — j3 from above. We have 
1 \ /_ 1 \ m 1 



1=1 



2 m _ m±i J y 2 m - / 2 m - V 2m+1 - (m + 1)/ 1 - |££ 

Since (1 + x) m > 1 + mx for all x € (-1, 1), it holds 

1 \ x 2 m + 1 -(m+l) _ 2 m + 1 -(m+l) T 2 m + ! _ £ ~ (m+1) 



l-fr^ [ 2 ™+i(l-p+l)] 2 



for all to > 3. Hence, q (l/(2 m — 2 2 1 )) > for all to > 3. If to = 2, the statement can be proved 
in the same way, just we use the exact expression (1 + x) 2 = 1 + 2x + x 2 instead of the estimate 
(1 + x) m > 1 + mx. 

Let us proceed to the extimate of 2 — [3 from below. 



1 \_( 2 1 V 1 | 1 

'■' I 2 m — m J y 2 m — — J 2™ — — — l — 



1 



2 m+1 -m ) V 2 m+1 
For all x G (-1, 0), it holds (1 + x) m < 1 + mx + (™)x 2 ; therefore, 

1 \ / to \ to m(rn — 1) m 

1-T-m - (l-T— rr) < 1-r—m + \, +77-I 1 



2 m+1 - my V 2" l+1 / 2 m+1 - to 2(2" l+1 - to) 2 2 m+1 

2 N 

111+2 -l 2to + to - 1 + 2 m+2 - 4to -l- — 



2(2 m +! - to) 2 V ' 2 

m 



1 - TO + — < 



2(2™+! - to) 2 V 2^ 
for all to > 2. Hence q (l/(2 m - f )) < 0. □ 
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Now we proceed to the eivenvalues /3j for j = 1, . . . , m — 1. For the sake of convenience let us 
set Bj := \(3j\ and jj := arg(/3j), i.e., 



Lemma A. 2. It holds 



= for all j = 1, . . . , to - 1 . 



|^.| <1 _ M5-4co S7j ) /_ln3 



for all j = 1, . . . ,m — 1. 

Proof. Since the value (3j = Bje ni is a solution of equation (|4"Tj) . necessarily 



(44) 



Hence 



|B™e im ^(2-B i e i ^)| =1 



S| m (4 - 4B 3 cos 73 + S?) = 1 



(45) 



Note that if to ^ 1, then obviously Bj rj 1. Therefore, equation (|45l) can be expressed approxi- 
mately as 

Bf 1 (4-4 cos7j + 1) k. 1 for to > 1 . 
Consequently, for to 3> 1 we have 



ln(5-4cos 7 -) 



Bi 



ln(5-4cos T -) 



2to 



2m 



■ ln(5— 4 cos 7j) 



With regard to this approximation, let us set 



1 _ ln(5 - 4cos7j) 



2m 



(1 + ^0, 



1 M5-4COS7,) (4g) 
2m 



(47) 



for all to, where 5j compensates the error of the approximation f|46[) . Comparing the statement (|44l) 
with the definition of <5j , we shall prove that 

Sj > for all j = 1, . . . , m — 1 . 

TO 

We proceed by contradiction. Let there be a j G {1, . . . , m — 1} such that <5j < — — . (Note 
that necessarily 6j > —1, because /3,-'s are of moduli less than one.) For all x > a > 1, it holds 



- 1 



< e ^ = (e a ) 



i+- 



Since ; = 1 — ^ for a; = 2m and a = (1 + dj) ln(5 — 4COS7-,), we have 



(l+jjj In(5-4co 3 Tjj 



< 



^ (1+ ^) ln( 5-4co S7j) y + --< 1+ ^)-<— = (5 _ 4cos7i) 



y 1 + d 3){ 1 +2 m -(l + ^)ln(5-4co S7j .) J 



Our assumption on <5j implies <5j < 0, therefore 
(1 <5j) ln(5 — 4 cos7j) 



< 



ln(5 — 4 cos 7j ) 



2m — (1 + (5j ) In ( 5 — 4 cos 7/) 2to — ln(5 — 4 cos 77) ' 
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hence 

i < (5 - 4cos 7 ^ 1+ ^ 1+ ^S^) . (48) 
At the same time we have from equation (|45[) 

i_ =4 - 4Bj cos 7j + S| = 5 - 4cos 7j + (1 - B i )(4cos7j - 2) + (1 - Bjf 

3 

> 5 - 4cos7j + (1 - Bj)(Acos lo - 2) . (49) 

Putting inequalities (|4*B"]) and (|4"9")l together, we get 

n A ln(5-4cos-yj) \ 

(5 - 4cos7 3 ) U+ d ^l i+ 2 m -in(5-4co 3 -, 3) ; > 5 _ 4 C os 7j + (1 - B i )(4cos7 i - 2) ; 

hence 

(5 - 4cos 7j )^ (1+ ^ — ^-""^ >! + (!- j?,) 4 "* 7 ' - 2 ■ 
J J 5-4cos 7j - 

This gives, with regard to equation J37J), 

e (^+(i+^) J-f^-Telf^) ) M5-4cos 73 ) _ x > ln(5-4cos7j) 5 4cos 7 j -2 

2m J 5 - 4 cos 7,- ' 

Since <5j < — ^ < — ln (5~4cos7j) by assumption, it holds 

g J + (l + fr) M5-4cos7i) 

therefore, the exponent in (|50[) is negative (or zero). Moreover, a simple analysis of the exponent, 
using the fact 6j > — 1, leads to the inequality 

(5, + (1 + fo) - ln (5- 4cos7 3 ) \ _ 4cog > _ for ^ ^ 

J J 2m - ln(5 - 4 cos 7,-)/ 

The convexity of the exponential function implies 

x , e h - 1 
- 1 < — - — 2: 


for all b < x < 0. Therefore, the left hand side of ([50)1 obeys 

e (^(Wj) 5^S^ ) ln(5-4co S73 ) _ x 

j. -ino/ ln(5-4cos7j) \ . , c . N 

< r - ^ ^ + 1 + ~m „ ; In 5 - 4cos7j) 

In 9 \ Jy 2m - ln(5 - 4cos 7i )/ v JJ 



( Si + (1 + 6* h ln( '" 4C °; 7j) A MS - 4cos 7j ) 
V 2m — ln(5 — 4 cos 7? ) / 



91n9 

Inequality (|50l) together with this estimate imply 

Sj + (1 + Sj ) n M5- 4COS7,) \ ln(5 _ > ln(5- 4cos7 3 ) (1 + 4cos 7j - - 2 



91n9 \ J J 2m - ln(5 - 4 cos 7,)/ J 2m 5-4C0S7,-' 

We divide both sides by ln(5 — 4cos7j), which is allowed due to 7 j ^ (recall that j3j ^ (0, +00) 
for all j = 1, . . . , m — 1); hence 

r , fl , r \ ln(5-4cos 7j ) 9 In 9 1 + 6, 4cos 7j - -2 

J 1 ^2m-ln(5-4cos 7j ) 8 2m 5-4cos 7j ' ^ ' 
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For all jj e K, In (5 - 4 cos 7,-) < In 9 and 

4cos7,-2 , 3 ,32 

^ = -H > -1 + - = — ; 

5 — 4 cos 7j 5 — 4 cos jj 9 3 

therefore, with regard to inequality (|51l) . 

In 9 91n9 1 + <5,- -2 3 to. 9.., _ , 

(5, + (1 + SA > • — = (1 + Si) . 

J y 3, 2m-\n9 8 2m 3 8m 1 3) 

Consequently, 

1 3 . . 

+ —)8 J > 



hence 



2m-ln9 8m/ J 2m-ln9 8m 



c 2m — In 9 8m 

5j > 



1 



2m-ln9 1 8m 
In 3 



This is a contradiction with the assumption Sj < — ^ , because 

-— < + ^ for all m> 2. 

m 1 + 75 — V" 5 + ir- 

2m — In 9 8m 



Lemma A. 3. The arguments of f3j satisfy 



□ 



/ 2ir ir 2ir ir \ . 

7 3 eh---j- + - (52) 

\ m om m om / 

/or a/Z j = 1, . . . , m — 1. 

Proof. Equation (l4"T1) has m + 1 solutions, namely 1, /3 and . . . , [3 m -\. Therefore, it suffices to 
show that every sector 



, 2ir ir 2ir ir . . 

B>0, 7 e j — ,j h - — )} for j = l,...,m- 1 

m om m bm 



contains exactly one solution of equation (|41[) . 
Let 

/3 = Bc i7 



be a solution of (|41l) . i.e., 

B m e im 7 ^ _ Bc i 7 ^ = j _ 

Hence 

m-7 = — arg (2 — _Be 17 ) + 2jn for a certain jgZ. (53) 

We can obviously assume j 6 {0,1, ...,m— 1} without loss of generality. Since the solutions 1 and 
/3 of equation (1531) are obtained for 7 = 0, and, therefore, for j = 0, we prove the statement in two 
steps: 1. We demonstrate that equation (151)1 has exactly one solution for every j = 1, . . . ,m— 1. 2. 
We show that the solution corresponding to j belongs to the sector Sj for every j = 1 , . . . , m — 1 . 
It holds 

2 - Be' 7 = 2 - Scos 7 - LB sin 7, 

hence 

— £?sin7 —sin 7 



tan(arg(2 - Be 17 )) = 



2 - 5 cos 7 ^ - cos 7 
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Furthermore, B < 1 implies 2 — B cos 7 > 0, hence 

arg(2-£e i7 ) e (-tt/2, tt/2) , 

i.e., we can write 

/ ;~\ —sin 7 

arg (2 — _Be 7 ) = arctan -5 

jy — cos 7 

To sum up, equation (|53p is equivalent to 



(54) 



7717 — arctan 



sin 7 



7? cos 7 



2j7T. 



(55) 



For every j — 1, ... , m — 1, the left hand side £(7) = H17 — arctan j_"" 7 — of equation ([55]) . 

s cos 

regarded as a function of 7 with a fixed -B < 1, is continuous and satisfies 

= L(0) < 2jir < 2m?r = L( 27r ) . 

Also, a simple calculation gives 



L (7) = m 



■ cos 7 — 1 



1 



2 • -| cos 7 + 1 



> m 



2. _ 1 

B 1 



> m - 1 > 0. 



Consequently, equation (|55p has indeed exactly one solution for every j = 1, . . . ,m — 1. The 
solution satisfies mf — 2jn £ (— tt/2, n/2). With regard to the numbering ((3]), we conclude that 

2j7T 7T 2jlT 7T 



' ' m 2m ' to 2to 

Now we improve this estimate in order to prove 7/ £ <Sj. Since 2/i3j > 2 for all j = 1, . . . , m — 1, 
we have 



■ sin 7j 



s 7j 



< 



Hi 



2 — cos 7j 



It is easy to show that 



hence 



sin 7 



2 — cos 7 

sin 7j 



< -= for all 7 € R , 
v3 



arctan 



b7 -cos7j 



< arctan 



1 _ 7T 

~ 6 



By substituting estimate ((56)) into equation ([55]). we obtain statement ([52 



(56) 
□ 
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